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ABSTRACT 



We have determined precise stellar parameters and lithium abundances in a sample of 117 stars with basic properties very similar 
to the Sun. This sample selection reduces biasing effects and systematic errors in the analysis. We estimate the ages of our sample 
stars mainly from isochrone fitting but also from measurements of rotation period and X-ray luminosity and test the connection 
between lithium abundance, age, and stellar parameters. We find strong evidence for increasing lithium depletion with age. Our sample 
includes 14 stars that are known to host planets and it does not support recent claims that planet-host stars have experienced more 
lithium depletion than stars without planets. We find the solar lithium abundance normal for a star of its age, mass, and metallicity. 
Furthermore, we analyze published data for 82 stars that were reported to support an enhanced lithium depletion in planet hosts. We 
show that those stars in fact follow an age trend very similar to that found with our sample and that the presence of giant planets is 
not related to low lithium abundances. Finally, we discuss the systematic biases that led to the incorrect conclusion of an enhanced 
lithium depletion in planet-host stars. 

Key words. Sun: abundances - stars: abundances - stars: planetary systems 



1. Introduction 

The lithium abundances of solar-like stars in the solar neigh- 
borhood spread over more than two orders of magnitude, which 
is much large r than the range o f abundances seen for other el- 
ements (e.g., Red dv et"ai1l2003l) . The Sun, in particular, has a 
very low lithium a bundance compare d to many nearby solar 
analogs (e.g.. lLambert & Reddvl 120041) . Furthermore, the pho- 
tospheric solar lithium abundance is about 160 times lower 
than that measured in meteorites (log eu,e = 1 .0 5 + 0.10 dexfl 
log 

Qj.met — 3.26 + 0.05 dex; both values are from lAsplund et al.l 
2009). This difference between the current solar and protosolar 
■ values is not predicted by stan dard stellar evolution models (e.g., 
' iD'Antona & MazziteUilll984l) . 

The wide range of observed lithium abundances in 
nearby solar-like stars is most likely due to a depen- 
dency between log eg a nd the star's age and mas s (e.g. , 
Montalban & S chatzma nl 120001: Charbo mier& Talonl 2005; 
Xiong & D eng 2009j iDo Nascimento etahl 120091) . Lithium is 



easily destroyed by proton capture reactions in stellar inte- 
riors. Thus, if lithium is transported between the chemically 
mixed outer convection zone and deeper lying regions with 
temperatures that are high enough for lithium desttuction, the 
photospheric abundance will decrease with time. Diffusion 
probably contributes to the lowering of the surface lithium abun- 
dance throughout the main-sequence stage. This would explain 
why the photospheric solar abundance is much smaller than 



1 We use the standard notation log e x = log + 12, where n x and n H 
are the the number densities of element X and hydrogen, respectively. 



Also, for 
log e F , 



metallicities we use the common abbreviation [Fe/H] = 



log 6°. 



the meteoritic one. We expect an enhanced lithium depletion in 
stars with larger convection zones on the main sequence as well 
as in stars with a higher degree of differential rotation between 
the radiative core and the convective envelope (see below). The 
reason is that lithium is only depleted as it moves to deeper 
and therefore hotter regions of a star, where the temperature is 
high enough (abou t 2.5 million K) for proton capture (see, e.g., 
|PinsonneaulJll997h . 

Recently, it has been suggested that the presence of plan- 
ets around a star could affect the evolu tion of the photospheric 
lithium abundance (e.g jBouvie r 2008). A long-lasting star-disk 
interaction during the star's pre-main sequence phase could slow 
down the host-star's rotation and therefore increase the degree 
of differential rotation between the star's core and envelope. 
Rotationally-driven mixing is then enhanced, thus destroying 
more lithium than in stars without planets because fast rotators 
evolve with little core-envelope decoupling. Planet migration af- 
fects the star's angular momentum, which might also have an 
impact on log eu- Finally, the inges tion of a planet can increase 
the surface lithium a bundance (e.g jMontalban & Rebolof 2002; 
llsraehan et al.ll200ll) . 

The possibility of a lithium-planet con nection is subjec t 
of ongoing di s cussio n s. Recent work by iGonzalez] ([2 008 ). 
Gonzalez et al.l d2010l) . ICastro et al.l ([2008), and Israelia n et al.l 
d2009l) sugg ests a possible log -plan et dependency, whereas 
iRvanl d2000t) and lLuck & Heiterl (^006) find that stars with plan- 
ets s how the same lithium di stribution as the comparison field 
stars. Tak eda et al.l (120071120 101) describe the stellar angular mo- 
mentum as the crucial factor that determines the lithium abun- 
dance of solar-like stars and find that slow rotators show an en- 
hanced lithium depletion. Planets could be the reason for a slow 
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rotation, but they were not able to draw firm conclusions due to 
the low number of planet hosts in their sample and the fact that 
their use of the star's projected rotational velocity, v sin i, instead 
of measured rotation periods introduces additional uncertainty, 
since the inclination angle i is unknown. 

From a practical point of view, an enhanced lithium deple- 
tion in planet-hosts would be greatly beneficial for the search 
for extrasolar planets, because all known methods for extrasolar 
planet detection (e.g., radial velocity, transits, or microlensing) 
are very time consuming. With an enhanced lithium depletion, 
however, one could preselect planet-host candidates just by 
obtaining the stars' chemical composition. 

In this paper, we derive lithium abundances and ages for a 
sample of solar-type stars to examine whether there is a correla- 
tion between lithium and age as well as a possible connection be- 
tween lithium and planets. We also examine lithi um abu n dance s 
and ages for the solar analog sample of llsraelian et al.l (|2009), 
who claim to have found evidence for an enhanced lithium de- 
pletion in planet-host stars. 



2. Observations & analysis 

Our sample consists of 11 7 solar-like s tars selected 
from the Hipparcos catalo g dPerrvman et alj 1 1997b as in 
iMelendez & RamfrezI d2007l) . They where observed at the 
McDonald (Robert G. Tull coude spectrograph on the 2.7m 
Harlan J. Smith telescope; RGT), Las Campanas (MIKE 
spectrograph on the 6.5m Magellan Clay telescope), and La 
Silla (HARPS spectrograph on the 3.6m ESO te lescope) obser- 
vatori es. Our few solar twins observed at Keck ( Melendez et al.l 
2006) are not discussed here since they are already included in 
the McDonald sample. 

The RGT and MIKE data (sp e ctra a s well as stellar pa- 
rameters) are f r om Ramirez et aT| (I2009L hereafter R09) and 
MelendezTtaD (120091 120101 hereafter M09), respectively; two 
stars (HIP10215 and HIP79672) are part of both samples. 
HARPS spectra for 12 more stars were obtained from the ESO 
archive, while spectra for 6 other stars w ere obtained from the 
S 4 N database dAllende Prieto et al J 1200411. One of the objects 
(HIP80337) occurs in both the HARPS and the S 4 N samples 
(we use the HARPS parameters, because they have the smaller 
uncertainties), so that we have 17 additional stars. The spectra 
for these stars were analyzed in an identical fashion as in R09 
(see below). All sub-samples contain one or more solar refer- 
ence objects for normalization: R09 used the light reflected from 
the asteroids Ceres and Vesta, M09 used Vesta, and for the stars 
added in this work, spectra of asteroid Ceres, Jupiter's moon 
Ganymedes, and the Moon were used. TableQ]shows the specifi- 
cations of each sub-sample, where the last two lines describe the 
additional data from this work. All spectra have a signal-to-noise 
ratio (S/N) greater than about 200, which allows the determina- 
tion of lithium abundances as low as solar. Note that three stars 
in our sample are also included in M09 (HIP79672) and R09 
(HIP14614 and HIP42438). For the further analysis, we use the 
parameters with the smaller uncertainties. 

The HARPS and S 4 N data were analyzed by first measur- 
ing Fe i and Fe n equivalent widths (EWs), where a line list of 
45 iron lines (34 Fe i and 11 Fe n lines) within the wavelength 
range from 4445 to 8294 A was used. The Fe lines have a broad 



Table 1. Specifications for the different sub-samples. 



sample 


instrument/ 


wavelength 


spectral 


number 


name 


telescope 


coverage 


resolution 


of 






(in A) 


R = AA./A 


stars 


R09 


RGT / McDonald 


3800-9125 


60,000 


63 


M09 


MIKE / Magellan 


3400-10000 


65,000 


42 


this 


RGT / McDonald, 


3800-9125 


45,000- 


18 


work 


HARPS / ESO 


3783-6865 


110,000 





2 The Spectroscopic Survey of Stars in the Solar Neighbourhood 
(S 4 N) data and more detailed information can be found at 
http ://hebe. as.utexas.edu/ s4n/ 



range of excitation potentials (from ~ 0.1 to ~ 4.6 eV). The 
line list adopted is from Scott et al. (in preparation, see also 
lAsplund et al1 l2009). To calculate the iron abundances ([Fe/H]) 
from the Fe i/Fe n lines , we used the spectrum synthesis program 
MOOG dSnedenll 1973b and A TLAS 9 model atmospheres with- 
out convective overshoot (e.g jKuruczll 19931) . We iteratively as- 
signed the stellar parameters effective temperature, surface grav- 
ity, and microturbulent velocity by forcing simultaneously Fe i 
excitation equilibrium and Fe i/ Fe n ionization equilib rium. For 
a more detailed description see iRamfrez et alJ J2009). With the 
method described above, we derived the following mean errors: 
o-(T eff ) = 40 K, cr(logg) = 0.06 dex, and cr([Fe/H]) = 0.025 
dex. 

Stellar ages and masses were determined primarily from the 
location of stars on the theoretical HR-diagram (T e ff vs. \ogg) 
as compared to theoretical predictions based on stellar evolution 
calcu lations (isochr ones). We produced a grid of Y 2 isochrones 
(e.g jYi et al.l l2001) with steps of 0.01 dex in metallicity around 
the solar value. For each star, we computed the age probabil- 
ity distribution of all isochrone points included within a 3-cr ra- 
dius from the observed stellar parameters. The adopted mean 
age and l-cr Gaussian-like upper and lo wer limits were derived 
from the age pr obability distribution (e.g.,|Lachaume et al. 1999; 
iReddv et alJl2003h . Although isochrone ages of unevolved stars 
are typically unreliable, the high precision of our stellar param- 
eters allowed us to infer reasonably accurate absolute isochrone 
ages, even for stars as young as ~ 3 Gyr; relative ages are 
naturally even better determined. For most stars younger than 
about 3 Gyr only upper limits could be determined. For these 
younger stars, we adopted ages ba s ed on measur ements of rota- 
tional per iods dGaidos et alfeoOOl lBarnesll2007l) and X-ray lu- 
minosity ( Agiieros et al. 2009) along with the rotation-age rela- 
tion by Barnes (2007) and the X-ray luminosity-age relation by 
iGuinan & Engld d2009l) . In a few cases of stars in the intermedi- 
ate age region (2 -4 Gyr), an average of isochrone and rotational 
ages was determined to improve the age estimate. For stars with 
ve ry unreliable isochrone ages we used the activity-based ages 
bv lMamaiek & Hillenbra nd (2008) and lSaffe et all d2005l) . Our 
adopted ages and the methods to obtain them are listed in Table 

13 

Using our stellar parameters as well as those in R09 and 
M09, we derived the lithium abundances for all 117 stars with 
line synthesis using MOOG. For this purpose we generated a line 
list from 6697 to 6717 A, i.e. 10 A around the lithium doublet 
at 6707.8 A. The whole wavelength range was synthesized with 
MOOG, where hyperfine splitting and the Li-doublet were taken 
into account. Knowing the basic stellar parameters, we were able 
to fit the lithium line including the effects of the projected rota- 
tional velocity v sin i and the microturbulent and macroturbulent 
velocities. We derived lithium abundances with a mean error of 
cr = 0.05 dex for stars in which the lithium doublet was detected. 
Our mean of all solar spectra is log eu = 1.03 ± 0.04 dex. 



P. Baumann et al.: Lithium depletion in solar-like stars 3 

Tabl e 2. Age and lithiu m abundance of solar twins in open clusters of near solar metallicity. Data are from the compilation by 
ISestito & Randichl (120051) . 



Cluster 


Age in Gyr 


log eu 


0"(log £L,) 


[Fe/H] 


Source 


IC2602 &IC2391 


0.030 


2.9 


0.1 


-0.05 


Randich et al. (2001) 


Pleiades 


0.07 


2.8 


0.1 


-0.03 


Soderblom et al. (1993) 


Blanco 1 


0.10 


2.9 


0.1 


+0.04 


Ford et al. (2005) 


M34 (NGC1039) 


0.25 


2.8 


0.1 


+0.07 


Jones etal. (1997) 


Coma Berenices 


0.6 


2.4 


0.15 


-0.05 


Fordetal. (2001) 


NGC762 


2.0 


2.1 


0.1 


+0.01 


Sestito et al. (2004) 


M67 


3.9 


1.2 


0.5 


+0.05 


Pasauini et al. (2008) 



Initially, we derived Li abundances assuming line formation 
in LTE (local thermal equilibrium), in ID, hydrostatic, plane par- 
allel ATLAS 9 model atmospheres. Abundance corrections were 
thereafter applied to account for departures from LTE in the for- 
mation of the Li resonance line . The non-LTE m odeling proce- 
dure is the same as described in iLind et alJ (120091) . but extended 
to cover also super-solar metallicities. For our sample stars, the 
abundance corrections range from -0.03 dex to +0.08 dex, de- 
pending on the lithium line strength and exact stellar parame- 
ters. In stars for which the equivalent width of the lithium line 
is below ~ 100 mA, over-ionization of neutral lithium results 
in positive abundances corrections, increasing with increasing 
metallicity and decreasing effective temperature. When the line 
starts to become saturated, increased photon losses pushes the 
statistical equilibrium in the opposite direction, i.e. into over- 
recombin ation, resul ting in negative corrections for some stars 
(see lLind et al.l2009l for more details). The non-LTE corrections 
are very small in comparison to the full range covered in lithium 
abundance, and hence the qualitative results of this study are the 
same for lithium abundances inferred in LTE and non-LTE. Note 
that t he NLTE corrections are computed using MARCS mod- 
els (|Gustafsson et al. 2008). Our NLTE corrected solar lithium 
abundance is 1 .07 + 0.04, in good agreement with the 3D-NLTE 
estimate bv lAsplund et alJ d2009h . 



Our adopted stellar parameters and derived LTE and non- 
LTE lithium abundances are given in Table [5] Fig. Q] shows the 
good agreement between the three observational sub-samples, 
which reduces errors introduced by inhomogeneous data pro- 
cessing and therefore makes the analysis more reliable. It also 
is a proof of the consistently good quality of the data. 



We have also considered the lithium abundances of so- 
lar twins from 8 open clusters as listed in Table [2] Data are 
from the compilation by ISestito & Randich] d2005l) as shown i n 
Table with updated data for M67 bv iPasauini et alJ d2008l) . 
The age for M67 is taken from VandenB erg et alJ d2007l) and 
Yada v et alJ (12 008). the lithium abundances for M67 stars are 
from lCastro et alJ (l2010h . The clusters IC2602 and IC2391 are 
listed as one here, because their parameters are basically the 
same. We only used open clusters around solar metallicity (0.0 + 
0. 1 dex) and with reliable data. The solar twins that we select in 
open clusters are stars of one solar mass by definition, i.e. they 
are main sequence stars with 1M based on their effective tem- 
perature. We take into account the increase of the solar effective 
temperature with increasing age in the selection of stars from 
open clusters. 




5400 



5600 



5800 
T etf (in K) 



6000 



Fig. 1. NLTE lithium abundance, logg, and metallicity plotted 
against effective temp era ture. R09 and M09 stand for data from 
Ramir ez et al . (2009])and iMelendez et all d2009i l2010h . respec- 
tively; TW is data re-analyzed in this work. 



3. Results 

3.1. Lithium and age 

Using our sample of solar-like stars a clear lithium-age trend is 
observed (Fig. |2). The dependency is as expected: older stars 
show more lithium depletion. The Spearman correlation coefH- 
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2 4 6 8 10 12 14 
Age in Gyr 



Fig. 2. log eu vs. age for stars from the three observational sub- 
samples. Down-arrows denote upper limits. The Sun is marked 
with 0. 



cient is r lol — -0.6 1 . This trend becomes clearer when we restrict 
the sample to solar twins, as in Fig. [3] We define solar twins as 
stars with [Fe/H] = 0.0 + 0. 1 and M = ( 1 .00 ± 0.04)M o The stars 
from the open clusters given in Table|2]fit the trend in Fig.[3]very 
well. This is very important, because the ages of these clusters 
are well determined and the fact that they lie in the midst of the 
lithium vs. age trend of the field solar twins suggests that the 
ages we derived for individual stars are reliable. The Spearman 
correlation coefficient for the solar twin log eu-age trend includ- 
ing the open cluster data is r tw j n = -0.75. Another interesting 
thing to point out here is the fact that the Sun (marked with in 
the figures) fits the trend very well. This leads to the conclusion 
that the Sun does not have a particularly low lithium abundance 
compared to stars of similar age, mass, and metallicity. 

Fig. [3] also compares our observationa l data with model pre- 
dictions from lCharbonnel & Talonl (120051) for different initial ro- 
tational velocities of the stars. These hydrodynamical models 
give predictions for the evolution of surface Li abundance in 
solar-type stars, accounting self-consistently for element segre- 
gation and transport of angular momentum by rotation, gravity 
waves, and meridional circulation. They agree not only with the 
general lithium depletion trend observed by us, but it could also 
explain the relatively large scatter as a result of differences in 
initial stellar rotational velocities. 

3.2. Lithium and planets 

In Fig. |4] we plot lithium abundance against age, this time 
for a sample of metal-rich solar analogs. As metal-rich solar 
analogs we define stars with [Fe/H] = 0.25 + 0.15 and M 
= (1 .08 ±0.08)M o . We use those criteria because our sub-sample 




2 4 6 

Age in Gyr 



Fig. 3. log en vs. age for solar twins from R09, M09, TW and 
from the solar twins in solar metallicity open clusters. Note the 
different scale compared to Fig. [2] The solid fines are the pre- 
dicted values from the models by ICharbonnel & Talonl (2005) 
for different initial rotational velocities. 



of planet-hosts is biased towards those higher metallicities and 
masses. In this case we make a distinction between stars that 
are known to host planets (filled symbols) and those for which 
planets have not yet been detected (open symbols). 

We used a two-dimensional Kolmogorov-Smirnov (KS) test 
to measure the probability that the samples of metal-rich solar 
analogs with and without planets in Fig. [4] belong to the same 
parent population. Using a Monte Carlo simulation, we took into 
account the errors in lithium abundance and age by choosing 
random, normally distributed values within each variable's 1- 
cr environment on the linear scale, which means that instead of 
log eu, we used 10 logeLi ~ 12 , that is The upper limits were ac- 
counted for by distributing the values uniformly between and 
the upper limit. 

We averaged the outcome of 1,000 KS tests. This gave a 
probability for our metal-rich solar analogs with planets and 
those without planets to be part of the same parent sample of 
64 + 15%; if we ignore the error bars and upper limits, this prob- 
ability goes up to 80%. This is very important for the further 
analysis of the data, because it tells us that there is no intrin- 
sic difference between the two sub-samples. It is highly unlikely 
that the planet-hosts and comparison stars are different regarding 
their surface lithium abundance. 

The age-lithium correlation coefficient for the solar twins 
is similar to that corresponding to the metal-rich solar analogs 
('twill = -0.75, ^analog = -0.71). However, the shapes of those 
trends are not identical. In the range from 3 to 6 Gyr, in particu- 
lar, it is clear that for a given age, metal-rich solar analogs have 
on average lower lithium abundances than solar twins (see also 
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2 4 6 8 

Age in Gyr 

Fig. 4. Same as Fig.[3]but for metal-rich solar analogs ([Fe/H] = 
0.25 + 0.15, M= 1.08 + 0.08 M ). 

Fig. |5](c)). This is independent of whether the star has a planet 
or not. The age-lithium trend in Sun-like stars is thus metal- 
licity dependent. This metallicity effect is predicted by stellar 
mode ls due to the deep er convection zone in more metal-rich 
stars (Castro et al. 2009) and has lately been confirmed (see, e.g. 
Ido Na sciment o et al.ll2010l Fig. 5). Note, however, that the mass 
ranges being compared are somewhat different, and that this will 
have an impact on the lithium evolution of those two samples. 

4. Discussion 

Recently, it was claimed that planet-harboring solar-type stars 
show an enhanced lithium depletion and that lithium surface 
abundances in Sun-like stars do not co rrelate with stellar ages 
dlsraelian et alj|2009t ISousa et al.|[20ToL hereafter 109 and S10, 
respectively). 

In Fig. [5] we plot age versus lithium abundance for all 82 
stars used in 109 along with the objects from this work (here- 
after B103). With the stellar parameters adopted by 109 we de- 
rived the ages for that sample using the same techniques as for 
our sample; the ages we derive for the 109 sample are given 
in Table [5] Fig. [5] shows the results for all stars panel (a), the 
solar twins (panel (b)), and the metal-rich solar analogs (panel 
(c)). We are using the same selection criteria for solar twins and 
metal-rich solar analogs as in Sect. [3] The agreement between 
the age-lithium relation found with our sample and that by 109 is 
excellent, in particular when looking at the solar twin plot. This 
shows t hat the stellar para meters used by 109 (which were de- 
rived bv lSousa et alJl2008l) are essentially on the same scale as 




4 6 8 
Age in Gyr 



3 For consistency, we used our LTE lithium abundance in this discus- 
sion because the 109 work does not take into account non-LTE correc- 
tions. 



Fig. 5. Comparison between the sample from Israelia rTet al.l 
d2009l 109) and our sample (B10). The solid line, identically 
drawn in each panel, is an arbitrary reference line to guide the 
eye to the different log eu levels in solar twins an metal-rich solar 
analogs. Note that for consistency we use the LTE log eu values 
here. 



ours and therefore the combination of both samples for this anal- 
ysis does not introduce systematic errors. In fact, for the 10 stars 
in common between our sample and 109 we find differences of 
3 + 20 K in T eff , 0.02 + 0.04 in log g, 0.003 + 0.023 in [Fe/H], and 
0.06 + 0.1 1 in log en (the latter for the only 3 stars with lithium 
doublet detection). 

In Fig. |5{a), ten stars with ages greater than 4 Gyr and 
higher than average log eu can be seen above the main locus. 
Taking a closer look at those "outliers", the most interesting fact 
is that they have a particularly low surface gravity (logg ^ 4.1) 
compared to the rest of stars. In Fig. [6] we show log eu vs log g. 
The surface lithium abundance on the low-log g side does not 
follow the main track, for which log eu decreases with lower 
surface gravity, which is essentially the age effect, given that 
all these stars have similar masses. To exclude the possibility 
of systematic errors in the analysis of the low \ogg outliers, 
we compared the parameters and in particular their ages with 
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Fig. 6. log eu vs. logg for the stars from our sample (B10) and 
109. 



various sources (see Table[3]l. Our derived ages for these outliers 
are in reasonably good agreement with respect to the the values 
given in the literature. Only two stars appear to be older than the 
ages given in the consulted references, but even that difference 
is only around 2 Gyr and therefore not big enough for these stars 
to cease being outliers. This leads us to the conclusion that the 
ages we derived for the 109 sam ple and the stellar parameters 
adopted by 109 (mostly derived bv lSousa e~a l. 2008) are correct 
and the high-lithium envelope in the lithium-age plane is most 
likely real. Thus, when restricted to a narrow range of T e fr 
around the solar value, logg ^4.1 stars have higher lithium 
abundances than less evolved stars of similar age. 

We have examined the results by S10, who claim that 
there is no correlation between lithium and age in the 109 
sample. The S10 sample is basic ally the sam e as in 109, but 
reduced to the 60 stars studied in ISousa et al.l (2008) because 
of the high homogeneity of the stellar parameters. We show 
their lithium-age trends in Fig. [7] The selection criteria we 
used for the solar twins and the metal-rich solar analogs are 
the same as in Fig. However, this time we are using the 
masses and ages determined by S10. Although their full sample 
appears to show no correlation (Fig. |7^), the solar twin sample 
shows a clear trend between lithium and age. Note that the 
one planet-host in this sample has a low lithium abundance 
because of its old age and not the fact that it hosts a planet. 
There is no clear correlation between lithium and age for the 
metal-rich solar analog sample in this case, but this could be 
due to uncertain ages. Since the solar twin sample includes only 
6 stars, we define another sample of "extended solar twins" 
with [Fe/H] = 0.0 + 0.1 and M = 1.00 + O.1OM . The resulting 
figure shows a very definite trend of log eu with age and only 
a single outlier appears. This outlier (HD215456), however, 
shows a relatively low logg of 4.10 (and an almost solar mass 
ofl.04M o ). 

We have also examined the lithium vs. effective temperature 
diagram presented by 109. As shown in Fig.[8](a), they found that 
almost all stars with a high lithium abundance (log eu >1.5 dex) 
have not been shown to be planet hosts yet, whereas planet hosts 
and objects where no planets have been found are distributed 
quite equally below that lithium abundance, although the high 
number of upper limits makes a direct comparison difficult. In 
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Fig. 7. log eu vs. age for the S10 sample. The selection criteria 
for the four panels are given in the text. Masses and ages adopted 
to make this figure are from S10. 



order to make a more robust comparison, we have restricted the 
comparison sample using the following criteria: we only con- 
sidered comparison objects within a 2<x range in [Fe/H], \ogg, 
and T e ff around planet hosts, where cr are the averag e values of 
the un certainties in the stellar parameters given by ISousa et al.l 
(2008). In this way, we make sure that all stars lie within the 
same region of parameter space and are not influenced by the 
age or metallicity effects we find. Note that we do not restrict 
the lithium range, only metallicity, surface gravity and effective 
temperature. Using this selection allows for a homogeneous and 
unbiased comparison. When we restrict the parameter range cov- 
ered by the comparison stars as described above, the lithium- 
planet connection disappears; as seen in Fig.[8](b), it is not pos- 
sible to conclude on stronger lithium depletion in planet hosts 
compared to other stars. We stress that this figure is plotted di- 
rectly from the 109 data without further manipulation or use of 
new parameters. 

Three systematic biases have led 109 and S10 to conclude 
that solar-type planet-hosts feature an enhanced lithium deple- 
tion and that there is no age dependence: 

1 . At [Fe/H] =2 0.0, the existing HARPS sample of solar analogs 
with planets are on average older and therefore more de- 
pleted in lithium than non-planet-hosts, but not because they 
have planets. 
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Fig. 8. Lithium abundance as a function of T e g in stars with and 
without detected planets from the 109 sample. Filled and open 
circles represent stars with and without detected planets, respec- 
tively. In the upper panel the original comparison made is shown, 
which is not appropriate because the sample being compared 
span different ranges in evolutionary phases and metallicities. 
A less biased comparison is shown in the bottom panel, where 
we only plot stars without detected planets with stellar param- 
eters (Teff, logg, [Fe/H]) within 2cr of the planet-hosting stars. 
When a proper comparison is made, i.e., using stars with similar 
fundamental parameters, lithium is not abnormally low in stars 
with detected giant planets. 

2. At higher [Fe/H], where most of the 109 planet-hosts con- 
centrate, there is a slightly different log eu vs. age trend such 
that, at a given age in the 3-6 Gyr range, metal-rich solar 
analogs are more lithium-poor compared to solar metallicity 
ones. This is true for both planet-hosts as well as single stars. 

3. 109 and S10 samples include a number of peculiarly high 
lithium abundance and relatively low log g (=; 4. 1) stars; only 
one of them is a planet host. 

The apparently lower lithium abundances of planet-hosts found 
by 109 can thus be fully explained by a combination of age and 
metallicity effects, not separately but through the age vs. lithium 
relation. 



5. Conclusions 

1. In stars of solar mass and solar metallicity, it is clear that 
older stars have experienced more surface lithium depletion. 
Both the overall lithium-age trend as well as the scatter that 
we observe in this sampl e of stars can be explai ned by the 
theoretical models by ICharbonnel & Talon! (120051) . 

2. Metal-rich ([Fe/H]~ 0.25) solar analogs (M~ l.O8M ) also 
exhibit a lithium-age trend, which is different from that seen 
in 1M , [Fe/H]= 0.0 stars. At any given age in the 3 to 6 Gyr 
range, the metal-rich solar analogs are more lithium-poor. 
This is true for both planet-hosts and single stars. 

3. For solar-like stars, the lithium vs age trends for planet-hosts 
and stars where no planets have been found are statistically 



identical. Thus, the presence of a planet does not influence 
the observed surface lithium abundance. 
4. A number of solar-like stars with unusually high lithium 
abundance for their age are present in the field. We note that 
all of them have relatively low logg 4.1. We intend to 
pursue further observational work to better understand this 
small group of relatively low surface gravity and peculiarly 
high lithium abundance. 

Acknowledgements. We thank G. Israelian for sending us the lithium abundance 
data from 109. 
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Table 3. Ages, masses, and lithium abundances for the outliers in Fig.|5ja). R' HK denotes ages derived from chromospheric 
activity, rot denotes ages derived from rotation periods. 
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unknown 


88427 


- 


5810 


57 


4.42 


0.07 


-0.160 


0.025 


0.85 


0.87 


-1.00 


0.97 


0.01 


5.7 


1.5 


iso 


R09 


no 


89162 


165357 


5835 


30 


4.32 


0.06 


0.07 


0.03 


1.35 


1.39 


0.09 


1.05 


0.01 


6.0 


0.6 


iso 


M09 


no 


89443 


238838 


5796 


73 


4.48 


0.12 


-0.020 


0.038 


1.06 


1.09 


-1.00 


1.01 


0.01 


4.4 


1.9 


iso 


R09 


no 


89650 


167060 


5855 


25 


4.48 


0.05 


0.02 


0.02 


1.03 


1.06 


0.07 


1.05 


0.01 


2.2 


1.0 


iso 


M09 


no 


91332 


171918 


5775 


25 


4.20 


0.05 


0.206 


0.03 


1.40 


1.46 


0.09 


1.12 


0.02 


6.4 


0.7 


iso 


M09+V05 


unknown 


96402 


184768 


5713 


49 


4.33 


0.032 


-0.029 


0.030 


0.64 


0.68 


-1.00 


0.97 


0.01 


8.7 


0.7 


iso 


TW+T07 


no 


96901 


186427 


5737 


28 


4.34 


0.04 


0.055 


0.016 


1.38 


1.42 


-1.00 


1.00 


0.01 


7.5 


0.5 


iso 


R09+V05+LH06+T07 


yes 


96895 


186408 


5808 


39 


4.33 


0.05 


0.097 


0.020 


0.81 


0.86 


0.06 


1.05 


0.01 


6.0 


0.6 


iso 


R09+V05+LH06 


unknown 


100963 


195034 


5802 


17 


4.45 


0.03 


0.008 


0.013 


1.72 


1.75 


0.05 


1.02 


0.01 


3.6 


0.8 


iso 


R09+T07+T09 


no 


100970 


195019 


5823 


40 


4.23 


0.026 


0.083 


0.025 


1.39 


1.43 


0.10 


1.06 


0.01 


7.2 


0.4 


iso 


TW+V05 


yes 


102152 


197027 


5737 


47 


4.35 


0.06 


-0.010 


0.022 


0.65 


0.69 


-1.00 


0.98 


0.01 


7.7 


0.9 


iso 


R09+M09 


no 


104504 


201422 


5836 


48 


4.50 


0.06 


-0.160 


0.022 


2.42 


2.42 


0.06 


1.00 


0.02 


3.0 


1.4 


iso 


R09 


no 


107350 


206860 


6015 


50 


4.48 


0.07 


-0.020 


0.019 


2.91 


2.87 


0.07 


1.09 


0.01 


0.2 


0.0 


rot 


R09+V05 


unknown 


108708 


209096 


5875 


51 


4.51 


0.07 


0.150 


0.024 


2.42 


2.44 


0.06 


1.10 


0.01 


1.3 


0.6 


iso 


R09 


no 


108996 


209562 


5838 


56 


4.50 


0.08 


0.060 


0.027 


2.40 


2.42 


0.07 


1.05 


0.01 


2.3 


1.1 


iso 


R09 


no 


109110 


209779 


5817 


60 


4.46 


0.033 


0.062 


0.030 


2.51 


2.52 


0.01 


1.04 


0.01 


0.7 


0.1 


rot 


TW+V05+T07+B07 


unknown 


109931 


- 


5739 


74 


4.29 


0.08 


0.040 


0.026 


0.95 


1.00 


-1.00 


1.00 


0.01 


8.2 


0.9 


iso 


R09 


no 


113357 


217014 


5803 


47 


4.38 


0.05 


0.221 


0.017 


1.36 


1.42 


0.07 


1.09 


0.01 


3.9 


0.7 


iso 


R09+V05+LH06 


yes 


115604 


- 


5821 


51 


4.43 


0.06 


0.140 


0.019 


0.85 


0.90 


-1.00 


1.07 


0.01 


3.1 


1.2 


iso 


R09 


no 


118159 


224448 


5905 


44 


4.55 


0.07 


-0.010 


0.022 


2.69 


2.68 


0.05 


1.06 


0.01 


0.8 


0.4 


iso 


R09 


no 



The abbreviations of the sources in Tabled are the following: 

R09, M 09 and TW are from iRamfrez et atl 1200911. IMelendez et all 42009H IMelendez et al] 120101) and this work as be fore. M 06 is IMelendez et al. 
j2006]) . V05 is IValenti & Fischer! 120050 08 i s ISousa et al l 12008). B0 7 Barnesl 120071) . L H06 [Luck & Heiteri f2006l) . PM0 8 iPorto de Mello et al 
(2008), T07 a nd T09 arelTakeda et al] <2007h and Takeda & Taiitsu ( 200% respectively, GOO islGaidos et aljfeOOOl) and MR07 is lMelendez & Ramirez 
J2007ri . S05 is lSaffe et al] JIpOl) , G09 is lGuinan & Engle ( 2009), and M08 lMamaieF& Hillenbrand ( 2008|). 
The -1 in cr(log e^) denotes upper limits. 



P. Baumann et al.: Lithium depletion in solar-like stars 
Table 5. Sample used in 109. Masses and ages are from this work. 



Star name or HIP 


HD 


Mass 


tr(m) 


Age 


O-(T) 


Star name or HIP 


HD 


Mass 


cr(m) 


Age 


o-(r) 


WASP 5 


- 


0.99 


0.06 


7.9 


3.3 


52409 


92788 


1.08 


0.01 


3.8 


1.0 


XO-1 


- 


1.01 


0.01 


2.3 


1.2 


53837 


95521 


0.98 


0.01 


3.4 


1.1 


1499 


1461 


1.07 


0.01 


4.5 


0.6 


54287 


96423 


1.01 


0.01 


7.2 


0.6 


1954 


2071 


0.97 


0.01 


4.6 


1.1 


54400 


96700 


0.97 


0.01 


6.6 


0.6 


2021 


2151 


1.12 


0.08 


6.7 


1.4 


97998 


97998 


0.90 


0.01 


1.9 


0.6 


5339 


4307 


1.01 


0.01 


9.0 


0.4 


60081 


107148 


1.12 


0.01 


3.0 


0.6 


6455 


8406 


0.98 


0.01 


3.0 


0.6 


60729 


108309 


1.05 


0.01 


7.5 


0.3 


8798 


11505 


0.93 


0.01 


8.6 


0.6 


62345 


111031 


1.10 


0.01 


3.4 


0.8 


9381 


12387 


0.91 


0.01 


9.1 


1.7 


64408 


114613 


1.20 


0.02 


5.9 


0.3 


9683 


12661 


1.10 


0.03 


4.5 


1.3 


64459 


114729 


0.97 


0.01 


9.7 


0.2 


12048 


16141 


1.09 


0.01 


6.9 


0.3 


64550 


114853 


0.92 


0.01 


7.4 


0.9 


12186 


16417 


1.12 


0.01 


6.7 


0.2 


65036 


115585 


1.13 


0.03 


5.3 


0.5 


14501 


19467 


0.94 


0.01 


10.0 


0.3 


71683 


128620 


1.17 


0.07 


4.1 


1.5 


15442 


20619 


0.94 


0.01 


3.9 


1.2 


74500 


134987 


1.10 


0.02 


5.4 


0.5 


15330 


20766 


0.94 


0.02 


3.4 


1.7 


78330 


143114 


0.88 


0.01 


9.9 


0.8 


15527 


20782 


0.98 


0.01 


7.3 


0.3 


78459 


143761 


0.98 


0.02 


6.1 


2.6 


16365 


21938 


0.86 


0.01 


10.8 


0.7 


79524 


145809 


0.96 


0.01 


10.3 


0.3 


19925 


27063 


1.01 


0.01 


4.2 


1.2 


79672 


146233 


1.03 


0.01 


3.3 


0.8 


20625 


28471 


0.97 


0.01 


7.7 


0.3 


83906 


154962 


1.22 


0.03 


4.7 


0.8 


20677 


28701 


0.89 


0.01 


9.5 


0.5 


160691 


160691 


1.14 


0.02 


4.8 


0.3 


23627 


32724 


0.97 


0.01 


9.2 


0.3 


95962 


183658 


1.01 


0.01 


5.3 


0.7 


22504 


34449 


1.02 


0.01 


1.5 


0.8 


96901 


186427 


1.02 


0.02 


5.0 


1.9 


25670 


36152 


1.05 


0.01 


2.6 


0.9 


97336 


187123 


1.07 


0.01 


3.5 


1.5 


26737 


37962 


0.94 


0.01 


5.2 


1.8 


97769 


188015 


1.10 


0.02 


1.8 


0.9 


27435 


38858 


0.95 


0.01 


3.3 


0.7 


98959 


189567 


0.92 


0.01 


8.4 


0.4 


30243 


44420 


1.11 


0.01 


3.5 


0.6 


98589 


189625 


1.09 


0.01 


2.5 


1.0 


30104 


44594 


1.08 


0.00 


4.1 


0.5 


102664 


198075 


0.99 


0.01 


2.3 


1.0 


30476 


45289 


0.97 


0.00 


8.8 


0.3 


104903 


202206 


1.09 


0.01 


1.4 


0.6 


34065 


53705 


0.97 


0.01 


6.8 


2.3 


106006 


204313 


1.06 


0.01 


4.6 


0.5 


36512 


59711A 


0.96 


0.01 


5.3 


1.0 


108468 


208704 


0.99 


0.01 


6.6 


0.3 


39417 


66428 


1.09 


0.02 


5.8 


1.0 


109821 


210918 


0.96 


0.01 


8.2 


0.4 


43726 


76151 


1.05 


0.01 


1.5 


0.5 


110109 


211415 


0.96 


0.01 


6.5 


1.2 


43686 


76700 


1.17 


0.07 


4.5 


1.2 


112414 


215456 


1.04 


0.01 


8.4 


0.4 


44713 


78429 


1.02 


0.01 


7.0 


0.5 


113357 


217014 


1.08 


0.02 


3.4 


1.6 


44890 


78538 


1.01 


0.01 


2.5 


1.1 




219542 


1.04 


0.02 


4.6 


1.5 


44860 


78558 


0.85 


0.01 


12.5 


0.7 


115577 


220507 


0.98 


0.01 


9.3 


0.5 


44896 


78612 


0.96 


0.01 


9.4 


0.3 


116250 


221420 


1.29 


0.06 


4.7 


0.7 


46007 


81110 


1.11 


0.01 


0.4 


0.1 


116852 


222480 


1.15 


0.03 


5.6 


0.8 


49728 


88084 


0.97 


0.01 


6.2 


0.8 


116906 


222582 


0.99 


0.01 


6.7 


0.8 


50534 


89454 


1.03 


0.01 


3.0 


1.1 


117320 


223171 


1.09 


0.01 


6.7 


0.3 


52369 


92719 


1.01 


0.01 


1.6 


0.9 


118123 


224393 


0.92 


0.01 


3.6 


1.0 



